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We study global asymptotic properties of a continuous time Leslie-Gower food chain model. We 
construct a Lyapunov function which enables us to establish global asymptotic stability of the unique 
coexisting equilibrium state. 
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In his papers [2l,'3l, P.H. Leslie introduced a predator-prey model where both interacting species are assumed to grow 
according to the logistic law. That is both species grow with a rate that is initially (for small population) proportional 
to the population and is limited by a carrying capacity. The novel feature of this model is that, while the carrying 
capacity for the prey is a positive constant, the carrying capacity of the predator's environment is proportional to 
the prey population. This idea leads to a model that is quite different from the Lotka-Volterra predator-prey model. 
Leslie's model stresses the fact that there are upper limits to the rates of increase of both prey, H , and predator, 
P, which are not recognised in the Lotka-Volterra model. These upper limits can be approached under favourable 
conditions: for the predator, when the number of prey per predator is large; for the prey, when the number of predators 
(and perhaps the number of prey also) is small. Furthermore, the Leslie-Gower model does not posses the "screw 
symmetry" that is inherent in the Lotka-Volterra model. This model was initially studied by Leslie and Gower Q, 
and then by Pielow 6]. In the case of continuous time, these considerations lead to the differential equations 6, p. 
91] 
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Here H{t) and P{t) are the prey and the predator populations respectively; r and q are the growth rates of the prey 
and the predator respectively; a is the attack rate; 1/rh is the carrying capacity of the prey environment, and 1/gc is 
the efficiency of consumption for the predator (that is H{t)/qc is the predator population that the prey population of 
the size H can support). All the constants in the system ([T]) are positive. This model always has the unique coexisting 
fixed point Q* ^{H*,P*), where 

^* = ^TT' = ^TTT' (2) 
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which was proved to be globally asymptotically stable 

The Leslie-Gower model can be immediately extended to the case of a food chain. The food chain composed of 
n+1 levels where the ith level depends (predates) upon only the i — 1 level can be represented by the transfer diagram 

H ^P,^P2^...^Pr.. 

This food chain can be described by the following system of differential equations: 

dH 
~dt 



{r -bH - aPi)H, 



f Pi PIP 



'^^^ = f 9i - Ci-^ SiPi+i \ Pi, (3) 



dt V"' 'P,_i 



dt V" "P„_ 



* First published in Mathematical Biosciences and Engineering 6:585—590 (2009). 
t andrei.korobeiniko v@ul.iel 
■t . William. lee@ul.ie. 



2 



Here the parameters qi, Ci and Si are defined by analogy to the single predator case, namely qi is the reproduction 
rate of the zth predator, q is defined so that 1/qiCi is the efficiency of consumption for the iih predator, and Si is the 
attack rate by the ith predator. In order for the equations to be biologically meaningful these parameters must all be 
positive quantities. 

The global properties of this model are given by the following Theorem: 

Theorem 1 . The Leslie- Gower chain model ^ always has a positive ( coexisting ) equilibrium state Q* = (if * , Pj* , . . . , P* ) ; 
this equilibrium state is unique and globally asymptotically stable. 

Proof. (1). Existence of the positive equiUbrium state. We prove this by induction. The positive equilibrium 
state always exists for n = 1: the coordinates of the equilibrium state are given by equalities ([2]). We assume that the 
statement of Theorem holds when n = m and prove that it holds for n — m + 1 as well. 

Starting from an m level chain in which, by assumption, all equilibrium populations ii(™)*, ^ i = 1, . . . , m are 

positive, we convert the system to an m + 1 level chain by introducing a population of top level {m + 1) predators and 
allow the system to equilibrate. It can readily be seen that the positive region IR.™+^ is an invariant set of this system. 
That prevents the sign of Pi for alH = 1, . . . , m changing. Thus, from the assumption of a positive equilibrium in the 
m level system, it follows that in the m + 1 level system all pj-"^~^^">* for z = 1, . . . ,m are positive or zero. Consider 
the final differential equation 

dPm+i , hPm+i , , 

In order for equilibrium population at P^+i = (rather than P^j^i > 0) to be attained we must have 

hp* 

lim a ^ < 0. (5) 

p^n - 

This requires limp^ ^_j.o P^^ < 0, whereas by assumption that the m- level system has a positive equilibrium state the 
converse holds. Thus the existence of an m level positive equilibrium state implies the existence of an m + 1 level 
positive equilibrium state. This completes this section of the proof. 

(2). Global asymptotic stability of the positive equilibrium state. A Lyapunov function 

ViH. P„ . . . , P„) = (in A + ^) + X: 5. (in ^ + f ) . 

where BiCi — Bi^iSi^iP*_^ and BiCi = aH*, is defined and continuous for all i7, Pi,...,P„ > 0. The function 
V{H, Pi , . . . , P„) satisfies 

dH H\ H )' dP, P\ P 

and hence the fixed point Q* is the only extremum of this function. It is easy to see that the point Q* is the global 
minimum of V{H, Pi, . . . , P„) in R"^^. (Fig. 1 shows the level curves of this function for n 
The function V{H, Pi, ... , P„) satisfies 

^=r-bH-aP,-r^+bH* + a^ 
dt H H 
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FIG. 1. Level curves of the Lyapunov function In + ^J- + B ^In -pi- + • 



By the definition of Bi, the equalities 



a— ^-BiCi— = 0, l^BiSi^ Z^^*c,-^=0 



hold. Furthermore, recollecting that 



P* P* 
r = bH*+aP^, ^ Ci—^ + s,P*+^, = Cn^^ 



hold at Q* , we obtain 



^ Sift = i^ici^ + E ^iCi;^ + E = «A* + 2 Bi-iSi_ii^*, 

i=l i=2 »-l 1=1 i=2 



4 



and 



p* p* p* 

= «A*^+E^»^»7J^13- 
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Therefore, 
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That is, for this model 7^ < strictly holds for all H.Pt, . . . ,P„ > 0, except the fixed point Q* where 7^ = 0. 



Therefore, by the Lyapunov asymptotic stability theorem [5^ , the fixed point Q* is globally asymptotically stable. 

(3). Uniqueness of the positive equilibrium state. At any equilibrium state, ^ = must hold. For this 
model, however, the fixed point Q* is the only point in M"^^ where ^ = holds. 

This completes the proof. □ 

Apart from the positive equilibrium state Q* where all n + 1 species coexist, this system also has n equilibrium 
states (where fc = 0,l,...,n — 1), which corresponds to the reduced m-species food chains 

H ^Pl^P2^...^Pk. 

Biologically these correspond to the case in which some external intervention has reduced the population of the ith 
species to zero (we have proved above that this system is uniformly persistent, and hence that can never occur for 
this model via the natural evolution of the system) leading to the extinction of the i + 1th level species that feeds on 
species i, and then to all higher levels of the food chain. For each of these equilibrium states, H, Pi, . . . , P^ > while 
Pk+i = . . . — Pn — 0- Thus, Qq corresponds to the predator-free case and has the coordinates Ho = r/b, Pi — ... — 
Pn = 0; Ql coincides with the equilibrium state ([2]) of the two-species model ([ij. The following Corollary immediately 
follows from the Theorem: 

Corollary 1. Apart from the positive equilibrium state Q*, the system has n non-negative equilibrium states Q\ 
(where A: = 0, 1, . . . , n — 1 Each of these equilibrium states is unstable in R>g, but globally asymptotically stable in 
the k- dimensional invariant subspace = {H, Pi, . . . , P/j > 0; Pk+i = . . . = P„ = 0}. 

In conclusion, we have to note that, apart from the mentioned n + I equilibrium states that are located in the 
nonnegative region IR.>q, the system has other n — 1 points with the coordinates that satisfy the equalities 

Pi P- 
r = bH + aPi, q., = + SiP^+i , (?„ 



-tz— 1 rn—1 

Indeed, it is readily seen that this system of algebraic equations is equivalent to a polynomial of the degree n and 
that this system has no complex solutions. However, the existence of these equilibria do not contradict the Theorem 
since these points are located outside of the non- negative region M^c which is the phase space of the system. The 
origin is an unstable equilibrium state of the system as well. 

These results demonstrate both the practicality and the usefulness of performing a stability analysis on non-trivial 
ecosystem models. We have shown, using a Leslie-Gower food chain model as an example, that it is possible to 
enumerate and characterize the stability properties of all the equilibrium states of the model. 
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